Relaxation vs decoherence: Spin and current dynamics in the anisotropic 
Kondo model at finite bias and magnetic field 
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Using a nonequilibrium renormalization group method we study the real-time evolution of spin 
and current in the anisotropic Kondo model (both antiferromagnetic and ferromagnetic) at finite 
magnetic field ho and bias voltage V. We derive analytic expressions for all times in the weak- 
coupling regime ma.x{V,ho,l/t} ^ Tc {Tc = strong coupling scale). We find that all observables 
decay both with the spin relaxation and decoherence rates ri/2. Various l/-dependent logarithmic, 
oscillatory, and power-law contributions are predicted. The low-energy cutoff of logarithmic terms is 
generically identified by the difference of transport decay rates. For small times t <C max{V, ho}~^ , 
we obtain universal dynamics for spin and current. 
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The real-time dynamics of small strongly interact- 
ing quantum systems coupled to several reservoirs (e.g. 
quantum dots, quantum impurities, or single molecules) 
is a fundamental nonequilibrium problem. The inter- 
est in this field stems from the experimental progress in 
controlling spin dynamics in quantum dot&l, and from 
the necessity to identify the qualitative dynamics for 
error correction schemes^. The theoretical description 
of such situations remains a huge challenge. Numeri- 
cal techniques like time-dependent numerical'^'^ and den- 
sity matrix renormalization group methods-, iterative 
path-integral methods^., and nonequilibrium Monte Carlo 
simulations! have been developed to describe the time- 
evolution. However, the description of finite bias or the 
long-time limit is often difficult. Analytically very little 
is known, except for special models which can be solved 
exactly^. For this reason, perturbative renormalization 
group (RG) methods have been developed for nonequi- 
librium problems^rJ^ to obtain results in the regime of 
weak coupling between dot and reservoirs. Concerning 
the time evolution, these methods have so far been ap- 
plied to the spin boson model^-ilfi and to the Kondo 
model at zero voltage for special regimesJii?. 

In this Letter, we will use a recently developed real- 
time RG method in frequency-space^iMii^ to calculate 
analytically the full time-evolution of the anisotropic 
Kondo model at finite voltage V and magnetic field Hq. 
This method has the particular advantage that the de- 
rived exact hierarchy of RG equations contains the full 
dynamics of local observables. In order to solve the RG 
equations systematically by expanding in the renormal- 
ized couplings, we will consider the weak-coupling regime 
where at least one physical low-energy scale is larger than 
the strong coupling scale Tc- For the Kondo model at zero 
temperature we consider 



Af = max{F,/io,lA} » T^, 



(1) 



where Tc = Tk (Kondo temperature) for the antiferro- 
magnetic (AFM) and Tc — for the ferromagnetic (FM) 
Kondo model. In this regime, the couplings J = J(At) 
can be chosen as expansion parameters, where J(A) de- 
notes the poor man scaling (PMS) solution at scale A 



(see ([3]) below). We find several interesting results which 
are proposed to be generic for any quantum dot in the 
Coulomb blockade regime where spin/orbital fluctuations 
dominate: 1. The voltage is an important energy scale 
for the dynamics which shows up in oscillatory, power- 
law, or logarithmic behavior. 2. In the long-time limit 
t ^ 1/ maxjy, /iq}, we find generically, i.e. in all orders 
of perturbation theory, that all terms are exponentially 
decaying with the transport rates F^. Furthermore, due 
to non-Markovian effects, each spin component and the 
current contain a sum of two terms, one decaying with 
the relaxation rate Fi and the other with the decoher- 
ence rate F2. 3. At finite voltage, resonances V,h^ 
\h — V\ -^OorV^h^O {h is the renormalized mag- 
netic field) are possible in the weak-coupling regime. In 
the limit l/max{y, ho} < i < l/\d\ (with S = h -V,h), 
we find logarithmic terms ~ (St) hi\dt\ for the transverse 
spin. In contrast to stationary quantitiesiiii^ii^ the cut- 
off scale of the logarithmic terms 6 is generically 
determined by the difference of decay rates Fi — F2. 4. 
In the short-time limit t <^ 1/ maxjV, /iq}, we obtain 
universal dynamics for spin and current. 

Model. — The model consists of a single spin-^, which 
is coupled by longitudinal and transverse exchange cou- 
plings J^^^ to the spins of two noninteracting reservoirs, 
see Fig. [l] for a sketch of the system. Experimental re- 
alizations of the model are provided by quantum dots^S 
and molecular magnets^!. The Hamiltonian reads H = 

Hres + Hd+ Hex, whcrC Hres = J2acrk <^aka\c,k°'aak dc- 

scribes two noninteracting reservoirs labeled hy a — L,R 
(cr =i,t denotes the spin and k is the state index), 
Hd — hoSz is the Hamiltonian of the local spin with 
Zeeman splitting, and 



He 



J1°' [SxSx + SySy) + jjO) 



(2) 



denotes the coupling^S. We use the notation s = 

5 J2aaka'a'k[ oLfc 3^a' o-a'a'k'^ (o. are the Paufi matriccs). 
The reservoirs are kept at different chemical potentials 
ML/-R = ±y/2 and are assumed to have a flat density of 
states in the band of width 2D. Furthermore, we consider 
the most interesting case of zero temperature. 




FIG. 1: (Color online) A spin-^ quantum system coupled via 
exchange couplings J^/x to the spins of two reservoirs. 



Method. — To study the dynamics of the local spin 
S.{t) = {S}{t) and the charge current I{t) = {i)it), we 
switch on the coupling 77ea; suddenly at the initial time 
t = 0. It means that we prepare the initial density ma- 
trix in the product form p(0) = P^j^ PlPRi where p^^ is 
an arbitrary density matrix for the local spin, and 
are grand canonical distributions for the left and right 
reservoirs. For t > 0, we calculate the dynamics from a 
kinetic equation for the reduced density matrix poit) — 
Trresp(t). Following Ref. H, one obtains in Laplace space 
and in Liouvillian notation pd{z) — i/{z ~ L'j^^ (z)) p^"*, 
where Vj^^ {z) is an effective dot Liouvillian. The spin 
dynamics follows from S_{t) = Trn S_pD{t), with poit) = 

J dz e~'^^* pd{z). The current in Laplace space is given 
by I{z) = —iTrDT,i{z)p]j{z), where denotes the 

current kernel. V^j^^ {z) and ^/(z) have been calculated 
in Ref. [13 for an arbitrary quantum dot in the Coulomb 
blockade regime with explicit formulas for the anisotropic 
Kondo model. The result is given by a systematic ex- 
pansion in the dot-reservoir couplings from PMS cut off 
at scale = max{V, /iq, 1^1} > Tc. The PMS equa- 
tions are well-known^^ and given by {d / dK) J ^^{K) — 
-2J_i(A) J^/2(A)/A, with the solution 



J. (A) = c 



1 



'Tj_ 
y A 



Ji(A) = 2c- 



1 



'i±. 

\ A 



(3) 



where = — J\ and Tr- are two invariants. In the 
isotropic case one obtains Jz/^{h) = 1/(2 ln(A/rft:)). 
For a study of the time dynamics it appears to be more 
convenient to expand in the couplings cut off at the time- 
dependent scale At defined in ([Ij. We have proven that 
this expansion is well-defined in the weak-coupling regime 
A( 3> Tc and equivalent to the expansion in J{Kz). 

In conventional perturbation theory one approximates 
L'j^^{z) w L'j^^ {z = 0) (Markov approximation) and con- 
siders only the terms up to second order in the bare 



exchange couplings (Born approximation). This 

approach fails to describe the current dynamics as well 
as yields simple exponential decay for the longitudinal 
(transverse) spin with the spin relaxation (decoherence) 
rate Fi (F2). In contrast, in this Letter we include non- 
Markovian terms and use a systematic perturbative ap- 
proach in the renormalized couplings Jz/_l- In partic- 
ular, replacing the effective Liouvillian in the resolvent 
l/(z — i^^(z)) by any of its nonzero eigenvalues Ap(z) 
evaluated up to O(J^), we obtain an expression of the 
generic form 



Ap + ^CbJb {z - Zb) In ■ 



A* 



1 -1 



-i{z - Zb) 



(4) 



The pole of this resolvent is denoted by Zp = Xp{zp). 
For the Kondo model, there are three nonzero poles 
{p ~ 1,±) at z\ = —iTi and z± = ±/i — iT2, where 
Fi = nJlih + nmx{V, h}) and F2 = {tt/2)J^V + Fi/2. 
The real parts of these poles up to 0{J) determine the 

quantities Ap, i.e. A± — ±(1 — Jz + ji^^)ho and Ai — 0. 
Besides the appearance of renormalized quantities, there 
are two new important contributions in (j?]). Firstly, 

the factor Z = I — 2{Jz — ji^"*) occurs, which contains 
Jz = Jz{At) and thereby determines the universal non- 
exponential short-time behavior. Secondly, the logarith- 
mic part is also of non-Markovian form. It induces branch 
cuts in the complex plane with branch points Zb located 
at Zb = ±h + nV — iT2 (for p — 1) and Zb = nV — iVi 
or Zb = ±/i + nV — iT2 (for p = ±), with 71 = 0, ±1 
(cf, and Jb = Jz/± are the certain coefficients and cou- 
plings). The branch cuts lead to exponential behavior 
on the scale z^, not Zp, i.e. with oscillation frequen- 
cies involving the voltage as well as unexpected decay 
rates. These exponential terms are multiplied by addi- 
tional power-law or logarithmic functions involving the 
scale Zb — Zp. Therefore, all logarithmic terms appearing 
in the long-time evolution are cut off by the difference of 
decay rates Tpp> = Fp — Fp' , a feature which is generic for 
all quantities entering the time evolution (cf . Ref. [13) . 

Results. — Considering all poles and branch cuts of the 
resolvents ^ and closing the integration contour of the 
inverse Laplace transform in the lower half plane of z, 
we obtain after some lengthy algebra the following final 
result for the spin and current dynamics in closed form 
(we set e = ?i = 1 and assume V,h>0) 



Sz{t) 

S+{t) 
lit) 



si'''' Ze-ri* + 5f (1 - e-^^') - Ji Re 



5: 



h+iri2 



j(0) ^-^.t^^ht 



Z - .Jl In(Ati) 



(0) 



IqV I i2 I Qh-iVx2 ^ Qh-V-iVi_-i 



JtsL^z^+v + Ji{s':-i[ - + ^r,rr-y") + -v) 



1 



-{Jl + 2Jl)V + TTji min{F, h}Sl\l - e~^'*) - - Si°^Jl Re 



j-h-v+iri2 



{V ^ -V) 



(5) 
(6) 
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FIG. 2: (Color online) Sz{t) in the isotropic Kondo model 
for ho = 100 Tk and various values of the applied voltage V, 
with Sz(0) = 1/2. Inset: S^(t) for V = ho ^ 100 Tk and 

^(O) = 1/2, Sy{0)=0. 
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FIG. 3: (Color online) Current I{t) in the isotropic Kondo 
model for ho = 100 Tk (solid lines) and ho = 200 Tk (dashed 
lines) and various values of the applied voltage V, with 
<S'z(0) — 1/2. We observe oscillations at times Ta-* ~ 0.05; 
the stationary current is reached for Tut > 0.5. 



where S+ = Sx + iSy and = Fi — The initial and 
the stationary magnetizations are denoted by 5^°^ = 5(0) 
and Sf = —h/{h + max{V, h}), and we have introduced 



the auxiliary functions 5"^, (t) — 
and If^, (t) = e 



In^+e 



-iz't 



Fa (at) 

iailn^ + e-*^'*Fi(at), with Fk{x) = 



Equations (IS|)-(I7]) provide the time dynamics up to 
O(J^) in all crossover regimes. We stress that all ob- 
servables decay both with the relaxation and decoher- 
ence rates. Correction terms ~ O(J^) oscillate with fre- 
quencies h±V , and, additionally, with h [nV] n = 0, ±) 
for the longitudinal (transverse) spin. From an exper- 
imental point of view these oscillations along with the 
occurrence of both decay rates are our most important 
results. Examples for the time evolution in the isotropic 
case Jz = J± are shown in Figs. [2] and [3l We now illus- 
trate the results in different time regimes. 

Long-time limit, off-resonance. — For large times t ^ 
max{V, /io}~^, we have At = Ac = max{y,/io} and 
Jz/± = Jz/±{^c)- In the off-resonance case V, h, \ h±V\ ^ 
Ac, we get typical power-law behavior from the asymp- 
totic expansions ^1(0;) « —i/x and F2(x) « 1/a;^, for 
a; I » 1. Besides the exponential decay, this gives addi- 
tional factors ~ l/{at)'', with fc = 1 for the current and 
k = 2 for the spin. For the transverse spin we obtain, for 
example. 



s+{t) 

-s 



J: 



.cos{Vt) 



(8) 



(0) Jl-Ttt 



JVt 



-iVt 



hH^ {h-V)H^ {h + V)H'^ 



with Z = Z - J'iHKt) + il{Jl + J2 [2 - 9(1/ - h)]]. 
To identify the different time scales from the exponential 
decay explicitly, situations where Fi and F2 differ sig- 
nificantly are of particular interest. For the anisotropic 
Kondo model at finite voltage this can be easily achieved 
since F2/F1 = 1/2 + 0{{Jz/Ji_f). Thus, for > 
J_L we obtain F2 » Fi. This is typically the case 
for single molecular magnets, where the transverse ex- 
change coupling is generated by quantum tunneling of 



magnetization^ (e.g. in Fe4 the ratio is J^f^ / J^^ ~ 
10~^). Fig. m shows the clear separation of time scales 
for Sx{t). The time scale at which the change of the 
leading behavior happens is approximately given by f ~ 
2 1 In J^l/Fai. At this time we also observe the change in 
the oscillation frequency from h to V . 

Long-time limit, on-resonance. — For large times t 3> 
A^^ and close to resonances, where |i5| <^ Ac, with 
5 — h ^ V 5 — h, we can enter the interesting time 
regime \5t\ <C 1. In this case, we obtain logarithmic con- 
tributions from Fk{at), with \at\ <C 1 and a = S — «Fi2. 
For the longitudinal spin and the current, the logarith- 
mic terms In |(5i| cancel out. In turn, for the transverse 
spin we obtain the logarithmic contributions 

S+ (<) ^ -a 6*^^/2 5(0) j2 g^ In I _ iTi2)t\, (9) 

with a = 1/2, 1 ior S = h — V, h. As a consequence, the 
spin perpendicular to the plane defined by the z— axis and 
the initial direction obtains a logarithmic enhancement 
at resonance S ^ after taking the derivative The 
cutoff of the logarithm at i5 = is determined by the 
difference F12 of the rates. Especially tor S = h V in 
the isotropic model the cutoff scale F12 = ^,Ph becomes 
zero close to resonance. In this case the perturbative 
approach breaks down for exponentially small values of 
ht where Jln|/it| ^ 0(1). This regime is beyond our 
approach and requires a nonperturbative treatment. 

Short-time limit. — For short timcs^"* t ^ A^^ 

we get At = 1/t and Jz/i_ = Jz/^{l/t) = J\i^^. The 
terms containing {S/ L)'^^, are negligible and we obtain 
from the first terms of ©-([Tl) the universal result 



[(j*)2 + 2(./i)2]. (11) 



S{t) = ZS!^"^ = (1-2J* + 2J(0))5(°\ (10) 
G{t) = ^ = ' 

For Ji°^ » and [(T^i)''''! < 1 (AFM case), we get 
m = S{0){1 + 2J(") - 2c - 4c{TKt)^'}, (12) 



(13) 
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FIG. 4: (Color online) logjQ(|5'i(f)|) in the anisotropic Kondo 
model for V = 2ho = 100 Tk and various values of the 
anisotropy = {ji°^f - {j'f^f, with S^O) = 1/2, S'y(O) = 
0. The dips have their origin in the oscillations of Sx{t). 

with c = \JiJi"^y - (ffy. For |(rj^t)4^| > 1 (FM 
case), we have to replace c — — c. As a result, univer- 
sal power laws are predicted in the short-time limit on 
the scale of the Kondo temperature. We note that for 
the FM Kondo model with V — ho — a power law 
has also been found for Sz{t) in Ref. HT). Here, we have 
found a complete analytic expression in terms of the PMS 
solution, together with results for the AFM model, the 
transverse spin, and the current. Furthermore, we iden- 
tified the validity range for finite by t ^ ^c^- 



isotropic case (c = 0), we obtain J* = — 1/(2 ln(Tft:i)) 

leading to S(t) = 5(0) (1 + 2ji°^ -I- l/ln{TKt)) and 
G{t) = (eV/i)37rV(81n2(TKt)) for both AFM and FM 
models. In this case, universal logarithmic terms occur, 
which have also been found for Sz{t) in Refs. l^flTl 

Summary. — We found complex dynamics for a spin 
coupled to fermionic reservoirs held at finite bias and for 
the current across it. The results go beyond Markovian 
theories and, in a well-controlled weak-coupling regime, 
are presented in closed analytic form covering all time 
regimes. Due to the generic form of Eq.(j3]), we propose 
our main conclusions to be generic for quantum dots in 
the Coulomb blockade regime. From Ref. ^141 it follows 
that the branch points of the logarithmic terms are gener- 
ically given by z^, ~ Zp + nV , where Zp is any pole of the 
leading order exponential decay, and n — 0, ±1, ±2, .... 
Therefore, we find in all orders of perturbation theory 
that all terms are exponentially decaying, all transport 
rates occur in principle in the decay of each observable, 
and the voltage appears in the oscillation frequencies. We 
have shown this explicitly for the Kondo model and re- 
vealed that the bias voltage is an important energy scale 
for the time evolution showing up in unexpected oscilla- 
tion frequencies and in various power-law and logarithmic 
contributions. 
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